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ABSTRACT 

Reduced  order  modeling  of 
thermoacoustic  instabilities  involves  the  coupled 
modeling  of  the  wave  propagation  in  the 
combustion  chamber  and  the  unsteady  heat 
release.  In  many  combustion  systems  the  Mach 
number  is  low  enough  that  the  effect  of  the  fluid 
motion  on  the  wave  propagation  can  be  ignored. 
Ignoring  the  fluid  motion  results  in  the  use  of  the 
wave  equation  to  model  the  wave  motion  in  the 
combustion  chamber.  In  a  previous  paper  the 
momentum  and  pressure  equations  were 
linearized  by  representing  the  fluid  motion  by  a 
steady  Mach  number.  In  that  research  the 
frequency  and  phase  relationship  change  as 
Mach  number  increases.  In  this  research 
unsteady  fluid  motion  is  considered.  The 
governing  equations  for  momentum  and 
pressure  are  modeled  in  SIMULINK  and  studied 
using  frequency  response  tools. 

1  INTRODUCTION 

Combustion  instability  manifests  itself  in 
many  of  today’s  combusting  system. 

Instabilities  have  been  reported  in  all  types  of 
propulsion  systems  including  rockets,  scramjets 
and  gas  turbine  combustors  and  augmentors. 
Throughout  the  50  year  development  of  high 
performance  gas  turbine  engines  combustion 
instability  has  always  been  a  major  challenge  for 
both  the  combustor  and  the  thrust  augmentor. 
Over  the  same  time  period  many  advances  in 
modeling  combustion  instability  have  been 
made. 

Traditional  thermoacoustics  modeling  is 
based  on  Culick’s  framework,  (Candel  2002, 
Malhotra  2004).  In  this  framework  pressure 
fluctuations  are  assumed  to  be  small  and  the 
Mach  number  of  the  flow  in  the  combustion 
chamber  are  assumed  to  be  small.  These 
assumptions  result  in  wave  equation  that 
describes  the  wave  propagation  in  the 


combustion  chamber  forced  by  an  unsteady 
heat  release  term.  This  type  of  model  has 
proven  very  effective  in  systems  where  the 
Mach  number  is  small. 

The  study  of  duct  flows  for  the  inlet  of 
gas  turbine  engines  has  shed  significant  light  on 
the  effect  of  steady  fluid  motion  on  wave 
propagation  in  ducts  where  fluid  motion  was 
appreciable,  Morse  and  Ingard  (1968), 

Eversman  (1970,  1971b),  Ingard  and  Singhai 
(1973,  1974),  Gogate  and  Munjal  (1992),  and 
Dokumaci  (1995).  The  conclusion  drawn  from 
much  of  this  research  is  the  wave  equation  does 
not  accurately  capture  the  wave  motion  in  the 
duct  when  there  is  appreciable  fluid  motion. 
Instead  a  more  complex  second  order  PDE  is 
required  to  capture  the  defect  of  the  fluid  motion. 

In  this  paper  the  steady  flow  the  effect  of 
unsteady  fluid  motion  is  studied.  The  wave 
propagation  is  model  using  finite  difference  in 
the  SIMULINK  environment.  The  steady 
pressure  and  momentum  equations  are  modeled 
using  finite  difference  to  reduce  them  to  a  set  of 
ODE’s.  These  ODEs  are  then  analyzed  using 
the  frequency  response  tools  available  in 
MATLAB.  Bode  plots  and  pole  zero  plots  are 
constructed  over  a  wide  range  of  operating 
conditions.  These  plots  are  then  compared  to 
previous  research  where  the  pressure  and 
momentum  equations  were  simplified  by 
assuming  the  fluid  motion  could  be  represented 
by  a  steady  Mach  number. 

2  WAVE  PROPAGATION  MODELS 

2.1.  Governing  Equations 

The  governing  equations  for  wave 
propagation  are  derived  from  the  continuity  and 
momentum  equations.  The  steady  flow  model  is 
derived  assuming  inviscid  flow,  with  no  body 
forces.  The  governing  equations  for  continuity 
and  momentum  reduce  to: 
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dj o  du  dp 

—  =  —p - u  — 

dt  dx  dx 

du  du  dp 

p —  +  pu —  =  — — 
dt  dx  dx 


(2-1) 

(2-2) 


Equation  2-1  is  further  reduced  by  assuming 
that,  for  an  isentropic  fluid,  pressure  and  density 
are  related  by: 

dp  =  a2  dp  (2-3) 


Equations  2-7  and  2-8  are  termed  the  Steady 
Flow  Model.  The  Homogenous  Wave  Model  is 
obtained  by  assuming  that  the  Mach  number  is 
zero.  Applying  this  assumption  results  in: 


P 


D  7  du 

-  -  -a  p — 
t  dx 

(2-9) 

du  dp 

dt  dx 

(2-10) 

The  pressure  equation  can  also  be  derived  from 
the  energy  equation  when  isentropic  flow  is 
assumed,  see  appendix  1.  Equation  2-1  then 
reduces  to: 


dp 

dt 


2  du  dp 

=  -a  p - a  — 

dx  dx 


(2-4) 


2.2  Wave  Propagation  Models 


2.3  Second  Order  Equations  for  Pressure 

The  equations  in  both  models  can  be 
combined  to  formulate  second  order  PDEs  with 
pressure  as  the  only  dependant  variable.  This  is 
accomplished  by  taking  the  temporal  derivative 
of  the  pressure  equation  and  the  spatial 
derivative  of  the  momentum  equation,  and  then 
substituting  the  derivative  momentum  equation 
into  the  derivative  pressure  equation.  Also 
required  is  the  assumptions: 


Equations  2-2  and  2-4  are  used  as  the 
governing  equations  for  both  the  Homogenous 
Wave  and  Steady  Flow  Models.  Equation  2-2 
can  also  be  modified  with  the  continuity  equation 
to  give  an  alternate  momentum  equation.  This 
methodology  was  applied  by  Morse  and  Ingard 
(1968)  to  achieve  a  second  order  PDE  that  was 
a  function  of  pressure  only.  See  Appendix  2  for 
the  derivation  of  Equation  2-5. 


du 

P^7  = 


u 


dt 


a 


dp 

dt 


+ 


uu 


\a 


-1 


^  dp 
j  dx 


(2-5) 


The  final  wave  propagation  models  are 
derived  by  assuming  the  flow  through  he  duct  is 
steady  and  thus  the  mean  is  constant.  The 
models  are  derived  by  assuming  the  velocity 
which  multiplies  the  derivatives  on  the  right  hand 
side  of  Equations  2-4  and  2-5  is  steady  and  can 
be  represented  with  Mach  number  by: 


u  =Ma 


(2-6) 


1 )  Constant  Properties 

2)  dp  »  dp 


For  the  Homogenous  Wave  Model  the  result  is: 


d2  p  2  d2u 

— P-  =  -a  p - 

dt2  dxdt 

d2u  _  d2 p 

^  dtdx  dx2 


(2-11) 

(2-12) 


Making  the  final  substitution: 


(2-13) 


Substituting  Equation  2-6  into  Equations  2-4  and 
2-5  results  in: 


dp 

dt 

P 


2  du 
=  -a  p  — 

dx 
du 
dt 


M  dp 
a  dt 


+ 


dx 


(2-7) 

(2-8) 


A  second  order  PDE  for  pressure  is  found  for 
the  Steady  Flow  Model,  Equations  2-7  and  2-8, 
by  applying  the  same  assumptions  and 
derivatives: 


d2  p  2  d2u  d2  p 

=  -a  p — --  Ma 


dt 


dx 


dx 


(2-14) 


2 


2007  AIAA  Aerospace  Sciences  Meeting  and  Exhibit 
8-11  January,  2007,  Reno  NV. 


AIAA  2007-0566 


dtdx  a  dtdx  ox1 


(2-15) 


Substituting  Equation  2-15  into  2-14,  the  second 
order  PDE  for  pressure  is: 


+  Cl2M 


d2p 

8x8t 


-  a 


(l  -  M 2 ) 


0 

(2-17) 


Equation  2-13  should  be  recognized  as  the 
Homogenous  Wave  Equation.  Equation  2-17  is 
often  termed  the  convective  wave  equation, 
Morse  and  Ingard  (1968)  and  Dowling  (2003) 

3.  SIMPLIFICATION  OF  PDEs  TO  ODE’s 
WITH  SECOND  ORDER  ACCURATE  FINITE 
DIFFERENCE 

The  first  order  equations  for  the 
Homogenous  Wave  Model  and  Steady  Flow 
Model  are  simplified  from  PDEs  to  ODEs  by 
using  finite  difference.  The  domain  is  broken 
down  into  n  nodes.  At  each  node  the  spatial 
derivatives  are  simplified  using  second  order 
accurate  finite  difference  operators.  For  the 
central  nodes  the  second  order  finite  difference 
operator  for  any  variable  v  is: 

dv  _v{x  +  Ax,y)-v(x-Ax,;y)  |  (3_i) 

dx  2Ax 


For  the  end  nodes  the  operators  are: 


dv 

dx 

dv 

dx 


-  3v(x,  y)  +  4v(x  +  Ax,  y )  -  v(x  +  2Ax,  y) 
2Ax 

3v(jc,  y)-  4v(x  -  A x,y)  +  v(x  -  2Ax ,  y) 
2Ax 


+  o[Ax2 ) 

+  o(Ax2 ) 


Applying  the  central  difference  to  the 
Homogenous  Wave  Model  Equations  2-9,  and 


2-10  become: 

dp  a2 p  /  \ 

sr~ 

(3-2) 

du  1  /  \ 

P^=  OA  'PM  Pi~  l) 
dt  2Ax 

(3-3) 

Applying  the  central  difference  to  the  Steady 
Flow  Model  results  in: 


dp 

dt 


“^  ((»)/+!  -  («),-!  )- 
2Ax 


Ma 

2Ax 


((p)m 


(pU  (3-4) 


ai  =  A/a,  +  k-i]((  , 

dt  a  dt  2  Ax 


(3-5) 


4.  BOUNDARY  CONDITIONS 

The  natural  frequencies  of  a  one 
dimensional  tube  will  differ  depending  on  the 
boundary  conditions  (Richardson  (1949),  Wood 
(1957)).  Depending  on  the  type  of  boundary  the 
pressure  and  velocity  will  either  be  a  node  or 
anti-node  at  the  boundary.  Since  acoustic 
pressure  and  acoustic  velocity  are  sinusoidal,  a 
node  is  defined  as  the  point  where  the  sinusoid 
is  zero.  An  anti-node  is  defined  as  the  point 
where  the  sinusoid  is  maximum.  Further,  if  a 
node  occurs  at  a  boundary  the  reflected  wave’s 
phase  will  be  preserved.  For  an  anti-node  at  a 
boundary,  the  reflected  wave’s  phase  changes 
by  180  degrees. 

For  ducts  without  flow,  there  are  two 
types  of  boundaries  in  the  extreme,  acoustically 
closed  and  acoustically  open.  For  an 
acoustically  closed  boundary  the  acoustic 
velocity  must  be  zero  at  the  boundary,  thus  the 
magnitude  of  the  acoustic  pressure  must  be 
maximum.  Since  the  waves  are  sinusoidal,  if 
the  acoustic  pressure  is  maximum,  the 
derivative  of  the  acoustic  pressure  must  be  zero 
(Pota,  and  Kelkar  (2001).  At  the  acoustically 
open  boundary  the  acoustic  pressure  must  be  a 
minimum  at  the  boundary,  thus  the  acoustic 
velocity  is  maximum.  Since  the  velocity  is  also 
sinusoidal,  if  the  acoustic  velocity  is  a  maximum, 
the  derivative  of  the  acoustic  velocity  must  be 
zero.  Applying  the  boundaries  to  the 
Homogenous  Wave  Equation  the  natural 
frequencies  can  be  obtained.  For  a  duct  with 
open-open  or  closed-closed  boundaries  the 
natural  frequencies  are: 


na 


For  a  duct  with  closed-open  or  open-closed 
boundary  the  natural  frequencies  are: 


W. 


c-o 

o-c 


(2 n  - 1  )a 
4L 


(4-2) 


5.  COMPUTATIONAL  DOMAIN 
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The  computational  domain  in  this  study 
is  a  one-dimensional  duct  whose  length  is  1.5 
meters.  At  the  inlet  of  the  duct,  on  the  left,  a 
closed  boundary  will  be  imposed.  At  the  exit  of 
the  duct,  on  the  right,  an  open  boundary  will  be 
imposed.  The  duct  will  contain  vitiated  air 
whose  temperature  is  approximately  1500K. 

This  temperature  results  in  a  wave  speed  of 
approximately  925  m/s  and  a  ratio  of  specific 
heats  of  1.3.  Table  5-1  is  a  list  of  the  first  10 
natural  frequencies  of  the  duct  assuming  no  fluid 
motion. 


N 

C-0 

Boundary 

1 

154.17 

2 

462.50 

3 

770.83 

4 

1079.17 

5 

1387.50 

6 

1695.83 

7 

2004.17 

8 

2312.50 

9 

2620.83 

10 

2929.17 

Table  5-1  Theoretical  Natural 
Frequencies  for  a  1.5  m  Closed-Open  Tube 
with  925  m/s  Wave  Speed 

6  COMPARISON  OF  HOMOGENOUS  WAVE 
AND  STEADY  FLOW  MODELS  FREQUENCY 
RESPONSE 


6.1  Homogenous  Wave  Model  Results 

The  computational  domain  was  modeled 
using  SIMULINK.  For  the  Homogenous  Wave 
Model  the  first  order  ODEs,  Equations  3-2  and 
3-3  were  applied  to  9  interior  nodes.  On  the  first 
and  last  nodes  the  boundary  conditions  were 
applied,  and  the  remaining  spatial  derivatives 
were  treated  with  the  upwind  and  downwind 
second  order  first  derivatives.  This  system  was 
then  analyzed  using  frequency-response 
methods  in  SIMULINK. 

Frequency-response  methods  are  the 
most  convenient  methods  for  engineers  to 
analyze  linear  systems.  In  these  methods,  the 
steady  state  response  of  a  sinusoidal  input  to 
a  system  is  found.  The  input  is  swept  over  a 
user  defined  range  of  frequencies.  The  results 
of  these  methods  are  reported  in  many  ways.  In 
this  research  the  results  of  the  frequency 


response  analysis  are  reported  with  semi- 
logarithmic  plots  called  Bode  magnitude  and 
phase  plots.  For  more  information  on  this 
method  and  Bode  plots  please  refer  to  Ogata 
(1970). 


Figure  6-1  is  a  plot  of  the  magnitude  of 
the  frequency  response,  or  Bode  Magnitude 
Plot,  of  the  Homogenous  Wave  Model  with 
closed-open  boundaries.  Unless  otherwise 
stated  the  input  point  is  always  the  first  node 
and  the  output  point  the  last  computational 
node.  Figure  6-2  is  the  Bode  Magnitude  plot  for 


M 


Figure  6-1  Bode  Magnitude  Plot  for  the 
Homogenous  Wave  Model,  Closed-Open 
Boundaries 


Bode  Magnitude  Plot  for  the  Homogenous  Wave  Model, 
Closed  -Closed  Boundaries 
3oo  11  Nodes 


Figure  6-2  Bode  Magnitude  Plot  for  the 
Homogenous  Wave  Model,  Closed-Closed 
Boundaries 

closed-closed  boundaries.  In  both  plots  there 
are  5  peaks  in  the  plot  at  5  distinct  frequencies. 
These  are  the  frequencies  where  the  model  is 
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most  sensitive  to  the  sinusoidal  input,  commonly 
referred  to  as  the  natural  frequencies. 

Table  6-1  lists  the  theoretical  and 
modeled  natural  frequencies  for  the 
Homogenous  Wave  Model.  It  lists  the  results  for 
both  a  closed-open  and  a  closed-closed  system. 
In  both  cases  the  finite  difference  model 
accurately  resolves  the  first  frequency  to  within 
less  than  1%.  For  the  higher  frequency  natural 
frequencies  the  mode  does  not  compare  a  well 
with  >15%  error  in  the  4th  natural  frequencies  for 
both  the  closed-open  and  closed-closed  models. 


1 

154.17 

154.10 

0.04 

154.1 

0.04 

2 

462.50 

448.00 

3.14 

461.7 

0.17 

3 

770.83 

705.00 

8.54 

767.2 

0.47 

4 

1079.17 

897.00 

16.88 

1069.2 

0.92 

5 

1387.50 

1365.5 

1.59 

6 

1695.83 

1657.5 

2.26 

7 

2004.17 

1940.1 

3.20 

8 

2312.50 

2215.8 

4.18 

9 

2620.83 

2480.5 

5.35 

10 

2929.17 

2733.9 

6.67 

n 

C-0 

11 

Nodes 

Error 

C-C 

11 

Nodes 

Error 

1 

154.17 

154.10 

0.04 

308.33 

308.00 

0.11 

2 

462.50 

448.00 

3.14 

616.67 

577.00 

6.43 

3 

770.83 

705.00 

8.54 

925.00 

795.00 

14.05 

4 

1079.17 

897.00 

16.88 

1233.33 

934.00 

24.27 

Table  6-1  Theoretical  and  Modeled  Natural 
Frequencies  for  1 1  Node  Closed-Open  and 
Closed-Closed  Homogenous  Wave  Models 

Figure  6-3  depicts  the  Bode  Magnitude 
Plot  for  the  Homogenous  Wave  Model  with  49 


Homogenous  Wave  Model,  Closed- 
Open  Boundaries,  49  Nodes 


Figure  6-3  Bode  Magnitude  Plot  for  the 
Homogenous  Wave  Model,  Closed-Open 
Boundaries,  49  Computational  Nodes 

computational  nodes.  Table  6-2  lists  the 
theoretical  frequencies  and  the  modeled  results 
for  1 1  and  49  computational  nodes.  Note  that 
as  the  number  of  nodes  increases  the  accuracy 
of  the  model  increases  tremendously.  For  the 
49  node  model  the  first  four  natural  frequencies 
are  predicted  to  within  less  than  1%  error. 


11 

49 

n 

C-0 

Nodes 

error 

Nodes 

error 

Table  6-2  Theoretical  and  Modeled  Natural 
Frequencies  for  Closed-Open  Homogenous 
Wave  Models 

6.2  Steady  Flow  Model  Results 

The  domain  was  modeled  employing  the 
Steady  Flow  Model  with  49  computational 
nodes.  In  this  model  the  fluid  motion  is 
considered  to  have  a  constant  mean  and  is  thus 
modeled  by  a  constant  wave  speed  and  Mach 
number.  Figure  6-4  depicts  the  Bode  magnitude 
plot  for  the  Homogenous  Wave  Model,  blue,  and 
The  Steady  Flow  Model  with  Mach  number 
equal  0.4,  purple.  Several  things  can  be  noted 
from  this  plot.  The  first  is  that  the  natural 
frequencies  are  shifted  to  the  left,  or  decrease 
by  over  1 5  Hz.  The  second  is  that  the 
magnitude  increases  over  40  decibels  (db). 
Finally  in  the  plot  for  M  =  0.4  there  is  a  second 
peak  with  lower  amplitude  a  few  hertz  larger 
than  the  first  for  each  natural  frequency. 

The  shift  of  the  natural  frequency  to  the 
right  is  addresses  first.  Figure  6-4  depicts  the 
Bode  Magnitude  plot  for  the  Homogenous  Wave 
Model  and  the  Steady  Flow  Model  with  Mach 
numbers  ranging  from  0.1  to  0.4.  When  the  first 
natural  frequency  at  each  Mach  number  is 
compared  to  the  Homogenous  Wave  Model  the 
frequencies  vary  by: 

t^nSF  ~  nHW  (l  —  M  ) 

This  result  can  be  explained  when  the 
result  is  considered  in  light  of  the  second  order 
PDEs  for  both  of  the  models,  repeated  here 
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dr 


dx1 


(2-13) 


~Y  +  a2M  — —  a2  (l  —  A/ 2 )— p-  =  0 

dr2  5x5?  v  ’  d2x 


(2-17) 


When  comparing  the  second  spatial  derivative  of 
each  equation,  the  wave  speed  in  the  Steady 
Flow  Model  is  modified  (1-M2).  Thus  one  would 
expect  that  the  natural  frequency  is  modified  by 
the  same  factor.  This  is  further  born  out  in 
Figure  6-5. 


Secondary  Peak  due  to  Finite  Difference 
Mach  Number  =  0.4 


Figure  6-4  Bode  Magnitude  Plot  for  the  Homogenous  Wave  Mode  (blue,  and  the  Steady  Flow 
Model,  (Purple)  with  Closed-Open  Boundaries,  49  Computational  Nodes 


In  figure  6-5  the  resultant  steady  state 
frequency  response  is  plotted  for  a  range  of 
Mach  numbers.  In  each  case  the  first  natural 
frequency  is  modified  by  (1-M2). 

The  secondary  peaks  in  Figures  6-4  and 
6-5  can  be  explained  when  finite  difference  is 
considered.  For  the  Homogenous  Wave  Model 
the  first  natural  frequency  is: 

a 

w  — - 

4 (n  -  lAx) 

For  the  Steady  Flow  Model  Figures  6-4  and  6-5 
indicate  the  new  natural  frequency  is 

Wn(mw}  -  4{n-l)Ax 


The  frequency  of  the  second  peak,  though,  is: 

ail-M2) 

W"2  ~  4(n  -  2)Ax 

The  second  peak  probably  occurs  because  a 
central  difference  is  used  on  the  center  9  points 
while  upwind/downwind  differences  are  used  on 
the  boundaries.  Thus  the  second  peak  is  the 
natural  frequency  if  the  central  nodes  are  only 
considered.  Another  explanation  for  the  second 
peak  could  be  a  reflection  from  one  of  the 
boundary  conditions. 

Figure  6-5  depicts  the  phase 
relationship  for  the  forward  (blue)  and  rearward 
(green)  waves.  The  “forward  traveling  wave” 
phase  relationship  is  obtained  by  using  the  inlet 
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as  the  input  and  outlet  as  the  output  for  the  frequency 


Figure  6-5  The  Effect  of  Mach  Number  on  the  Phase  of  the  Forward  and  Backward  Wave 

Propagation 

response  calculation.  The  green  line,  or  using  the  inlet  as  the  output  and  outlet  as 

“rearward  traveling  wave”,  is  obtained  by 
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the  input  for  the  frequency  response 
calculation. 

The  first  phase  plot  depicts  the 
Homogenous  Wave  Model,  M  =  0.0.  Note 
there  is  only  a  blue  curve  for  this  plot.  This 
is  because  the  forward  and  rearward  waves 
each  have  the  same  amplitude  and  phase 
relationship.  Mathematically  the  forward 
and  rearward  waves  combine  to  form  a 
standing  wave.  In  this  plot  the  phase  is  -180 
degrees  for  frequencies  less  than  the  natural 
frequency.  The  phase  then  increases  by 
1 80  degrees  at  the  natural  frequency.  It  is 
then  constant  again  until  the  next  natural 
frequency  where  it  again  increases  by  180. 

Equation  6-1  is  a  simple  second 
order  transfer  function.  The  frequency 
response  of 


this  system  is  displayed  in  the  Bode 
Magnitude/Phase  plot,  Figure  6-6.  For  this 
system  there  is  a  step  change  in  the  phase 
of  180  degrees  at  the  natural  frequency  of 
the  system.  The  phase  relationship  for  the 
Homogenous  Wave  Model  at  the  natural 
frequency  is  indicative  of  that  of  a  second 
order  system. 


Bode  Diagram 

From:  Ini  (pt.  1)  To:  Transfer  Fen  (pt.  1) 


10"2  10"1  10° 
Frequency  (Hz) 

Figure  6-6  Bode  Magnitude/Phase  Plot 
for  the  Simple  Second  Order  System, 

Eqn.  6-1 


The  second,  third  and  fourth  phase 
plots  are  for  Mach  numbers  0.1 , 0.25,  and 
0.4  respectively.  In  each  case  the  forward 
traveling  and  rearward  traveling  waves  have 
very  different  phase  relationships  than  the 
Homogenous  Model.  In  each  case  the 
phase  of  rearward  traveling  wave  is 
increased  by  180  degrees  relative  to  the 
forward  traveling  wave  at  low  frequency.  As 
the  input  frequency  is  increased  the  phase 
of  the  forward  traveling  wave  increases 
while  that  of  the  rearward  traveling  wave 
decreases.  Near  the  natural  frequency  the 
phase  of  the  forward  traveling  wave 
continues  to  increase  and  the  phase  of  the 
rearward  traveling  wave  changes  sign  and 
also  begins  to  increase.  At  the  natural 
frequency,  the  forward  traveling  wave’s 
phased  changes  by  -180  degrees,  while  the 
rearward  traveling  wave’s  frequency 
changes  by  -360  degrees. 

Equation  6-2  is  a  damped  second 
order  system  with  a  zero  and  two  poles  at 
the  origin. 

1  s  +  a 

—  ^  (6-2) 
s  s  +  lQs  +  a 

The  frequency  response  of  this  transfer 
function  is  depicted  in  Figure  6-7.  Figure  6- 
8  is  the  frequency  response  of  the  Steady 
Flow  Model  with  Mach  number  =  0.4.  Note 
the  phase  portion  of  the  forward  traveling 
wave,  blue,  Figure  6-8  has  the  same  shape 
as  the  phase  of  Figure  6-2 


Figure  6-7  Bode  Plot  of  Equation  6-2 

Equation  6-3  has  two  damped 
second  order  systems  with  similar  natural 
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frequency  multiplied  by  a  pole  and  zero. 
Figure  6-9  is 


s  +  b  s2  +  2£s  +  c  s2  +  2£s  +  d 
c  ~  d  >  a  >  b 

the  Bode  plot  for  Equation  6-3.  Comparing 
the  phase  plot  of  Figure  6-9  to  that  of  the 
backward  traveling  wave  in  Figure  6-8  they 
both  have  similar  shape.  In  both  cases  the 
phase  is  near  zero  for  low  frequency, 
decreases,  then  increases  as  the  frequency 
increases  to  the  natural  frequency.  At  the 
natural  frequency  the  phase  goes  through  - 
360. 


Figure  6-8  Bode  Plot  of  the  Steady  Flow 
Model  with  Mach  Number  0.4 


Figure  6-9  Bode  Plot  of  Equation  6-3 

From  this  analysis  three  things  can  be 
concluded  about  increasing  the  Mach 
number  in  the  Steady  Flow  Model 

1)  Increasing  Mach  number  has 
the  same  effect  on  the  phase  and  amplitude 
as  adding  damping  to  the  system. 


2)  For  the  forward  traveling  wave, 
to  match  the  behavior,  a  second  order 
system  near  the  origin  and  a  zero  are 
required  to  match  the  shape  of  the 
amplitude  and  phase  curves  of  the  Steady 
Flow  Model. 

3)  For  the  backward  traveling 
wave,  a  damped  second  order  system  near 
the  natural  frequency,  a  zero  and  a  pole  are 
required  to  match  the  shape  of  the 
amplitude  and  phase  curves  of  the  Steady 
Flow  Model 

7  CONCLUSIONS  AND 
RECOMMENDATIONS 

Steady  fluid  motion  has  a  three 
fold  effect  on  wave  propagation.  First, 
the  steady  fluid  motion,  represented  by 
Mach  number,  decreasing  the  natural 
frequency  of  the  system.  Second  the 
fluid  motion  changes  the  phase  behavior 
of  the  waves.  For  the  Homogenous  Wave 
Model  the  phase  is  -180  before  the 
natural  frequency  and  0  after,  with  a  180 
degree  phase  change  at  the  natural 
frequency.  For  the  Steady  Flow  model 
the  forward  and  backward  traveling 
waves  exhibit  variable  phase  through  the 
entire  range  of  frequencies.  Third  the 
amplitude  of  the  backward  wave  is  higher 
than  that  of  the  forward  wave. 

In  three  dimensional  cylindrical 
combustors  transverse  waves  will  tend 
to  spin  in  one  direction  instead  of 
remaining  steady.  For  cylindrical 
combustion  systems  with  swirl,  the 
different  phase  relationship  of  the 
forward  and  backward  waves,  as  well  as 
their  differing  amplitude  may  cause  the 
waves  to  prefer  to  spin  in  one  direction. 
According  to  the  modeled  data,  this 
spinning  behavior  may  even  be  true  for 
low  Mach  numbers,  ~0.1.  The  mode 
shapes  in  a  two  dimensional  chamber 
with  axial  and  transverse  flow  needs  to 
be  studied  further  to  determine  if  this 
hypothesis  is  true.  This  phenomena 
should  also  be  studied  in  the  presence  of 
loss  mechanisms  such  as  viscosity  and 
damping  added  to  the  model. 

A  second  peak  appears  when  the 
boundary  conditions  are  applied  to  the 
Steady  Flow  Model.  To  reproduce  the 
forced  response  with  transfer  functions 
second  order  systems  were  required  to 
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be  added  at  the  origin  or  near  the  natural 
frequency  of  the  system.  More  study  is 
required  to  determine  if  these  second 
order  systems  are  an  artifact  of  the  finite 
difference  simplification,  or  if  they  are  an 
artifact  of  the  boundary  conditions 
applied. 
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Appendix  1  Alternate  Pressure  Equation 

Starting  with  the  conservation  form 
of  the  internal  energy  equation: 

d{ep )  d(pue)  du 

^r+^r=-p&(A1-1) 

For  isentropic  flow  pressure,  density  and 
internal  energy  are  related  by: 

P  _  P  1  _ _ P 


P  = 


RT 


O'— i)’ 


pe  = 


r- 1 


Reducing  equation  A1-1  to: 

r  -  \  r  -  ^ 


d 


P 


\r- 


i 


d\ 


u 

v  Y 


P 


1 


du 

=  -p — (A1-2) 
dt  dx  dx 

Assuming  constant  properties, 
equation  A1-2  becomes: 

1  dp  1  d(up )  du 

y  - 1  dt  (y  - 1)  dx  dx 


+ 


Further  manipulations: 

d(up)  /.  \  du 

=  ~r)p— 

ox  ox 


dp 
~dt 

Applying  the  chain  rule  to  the  spatial 
derivative: 

dp  du  dp  \  du 

—  +  p  —  +  u  —  =  {\-  y)p  — 

dt  dx  dx  dx 

dp  dp  /  du 

—  +  u  —  =  {\-y-\)p  — 

dt  dx  dx 

The  final  pressure  equation  is: 

dp  dp  du 

—  +  u  —  =  -yp  — 

dt  dx  dx 


dp  dp 
—  +  u  — 
dt  dx 


-  fXl 


du 

dx 


(A1-3) 


du  du  dp 

p - F  pu —  = - 

dt  dx  dx 

Multiplying  equation  2-1  by  u  and 
rearranging  results  in: 


(2-2) 


pu 


du 


—  -u 


dp  ....  dp 


uu 


(A2-1 ) 


dx  dt  dx 

The  derivates  of  density  are  converted  by 
the  isentropic  relationship  for  pressure  and 
density: 

.2 ' 


dp  =  a  dp 

Resulting  in 

du 

pu-  =  -- 

Ox 


u  dp 
~~dt 


a 


uu  dp 
dx 


a 


(A2-2) 


The  alternate  momentum  equation  is 
obtained  by  substituting  Equation  A2-2  into 
Equation  2-2 

du  u  dp  uu  dp  dp 

dt  a2  dt  a2  dx  dx 

Grouping  terms  results  in: 

du  u  dp  ( ^ 

Por  =  ^—  + 1 


dt 


u  dp  (uu  A  dp 
a 2  dt  \a2  J  dx 


(A2-3) 


Appendix  2  Alternate  Momentum 
Equation 

Starting  with  the  momentum  and 
continuity  equations: 

dp  du  dp 

=  -p  —  ~U—  (2-1) 


dt 


—  u 

dx  dx 
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